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Non-equilibrium features of a first order phase transition from the quark- 
' gluon plasma to a hadronic gas in relativistic heavy-ion collisions are dis- 



cussed. It is demonstrated that strong collective expansion may lead to the 
fragmentation of the plasma phase into droplets surrounded by undersatu- 



< 

' rated hadronic gas. Subsequent hadronization of droplets will generate strong 

non-statistical fluctuations in the hadron rapidity distribution in individual 
events. The strongest fluctuations are expected in the vicinity of the phase 



> 

o 

C*~) ' transition threshold. 
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The main goal of present and future experiments with relativistic heavy ions is to produce 
and study in the laboratory a new form of strongly interacting matter, the Quark-Gluon 
Plasma (QGP). Due to the confinement of color charges, only colorless hadronic final states 



^ ■ can be observed experimentally. Therefore, QGP properties can be studied only indirectly 
through the final hadron distributions or by penetrating electromagnetic probes. Many 
QGP signatures have been proposed in the recent years, particularly ones which assume an 
equilibrium phase transition between QGP and hadronic gas. 

The phase structure of QCD is not yet fully understood. Reliable lattice calculations exist 
only for baryon-free matter where they predict a second order phase transition or crossover 
at T c « 160 MeV. Recent calculations using different models reveal the possibility of a 
first order phase transition at large baryon chemical potentials and moderate temperatures. 
The predicted phase diagram in the (T, fx) plane contains a first order transition line (below 
called the critical line) terminated at T fa 120 MeV by a (tri)critical point 0,|3|]. Possible 
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signatures of this point in heavy-ion collisions are discussed in Ref. ||. Under certain non- 
equilibrium conditions, a first order transition is also predicted for baryon-free matter ||. 

A striking feature of relativistic heavy-ion collisions, confirmed in many experiments (see 
e.g. [0]), is a very strong collective expansion of matter. The applicability of equilibrium 
concepts for describing phase transitions under such conditions becomes questionable. The 
goal of this paper is to demonstrate that non-equilibrium phase transitions in rapidly ex- 
panding matter can lead to interesting phenomena which, in a certain sense, can be even 
easier to observe. 

To make the discussion below more concrete, I adopt a picture of the chiral phase tran- 
sition for which the mean chiral field $ = (cr, tt) serves as an order parameter. It is assumed 
that the theory respects chiral symmetry, which is spontaneously broken in the vacuum 
where a = f n , ti = 0. The effective thermodynamic potential fl(T, /i; $) depends, besides 
$, on temperature T and baryon chemical potential fi. Since explicit symmetry breaking 
terms are supposed to be small, to a good approximation Q is a function of <£> 2 = a 2 + ti 2 . 

The schematic behaviour of f2(T, /i; <3>) as a function of the order parameter field o at 
7i = is shown in Fig. 1. The curves from bottom to top correspond to different stages of the 
isentropic expansion of homogeneous matter. Each curve represents a certain point on the 
(T, /i) trajectory. The minima of Q correspond to the stable or metastable states of matter 
under the condition of thermodynamical equilibrium, where the pressure is P = —Q min . 
The figure is based on the calculations within the linear a-model with constituent quarks 
|U, which predicts a rather weak first order phase transition. A similar structure of fi(T, fi; $) 
but, possibly, with a stronger phase transition is predicted by the NJL model 0] and by the 
random matrix model 0. The discussion below is quite general. 

Assume that at some early stage of the reaction thermal (but not necessarily chemical) 
equilibrium is established and partonic matter is in a "high energy density" phase Q. This 
state corresponds to the absolute minimum of Q with the order parameter close to zero, 
a rj 0, 7i ~ 0, and chiral symmetry restored (curve 1). Due to a very high internal pressure, 
Q matter will expand and cool down. At some stage a metastable minimum appears in 
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£7 at a finite value of a corresponding to a "low energy density" phase H, in which chiral 
symmetry is spontaneously broken. At some later time, the critical line in the (T, //) plane 
is crossed where the Q and H minima have equal depths, i.e. Ph = Pq (curve 2). At later 
times the H phase becomes more favourable (curve 3), but the two phases are still separated 
by the potential barrier. If the expansion of the Q phase continues until the barrier vanishes 
(curve 4), the system will find itself in an absolutely unstable state at a maximum of the 
thermodynamic potential. Therefore, it will freely roll down into the lower energy state 
corresponding to the H phase. This situation is known as a spinodal instability. 

According to the standard theory of homogeneous nucleation || , supercritical bubbles 
of the H phase appear only below the critical line. Under condition of thermal equilibrium 
between the two phases, the supercritical bubbles can only grow through the conversion of 
new portions of the Q matter into the H phase on the bubble boundary. The bubble growth 
is then limited by a small viscosity of the Q phase resulting in a slow dissipation of the 
latent heat ||. Therefore, a certain degree of supercooling is needed in order to convert a 
significant fraction of the Q matter into the H phase in the form of nucleation bubbles [|5],|5] . 

In rapidly expanding matter the nucleation picture might be very different. Let us 
consider first an isotropically expanding system with the collective velocity field which follows 
the Hubble law locally, v(r) = H -r. The Hubble "constant" H may in general be a function 
of time, e.g. H ~ l/t. Obviously, the same velocity field is seen from any local rest 
frame comoving with the matter. Suppose that a bubble of the H phase has formed in 
the expanding Q matter because of a statistical fluctuation. The change in thermodynamic 
potential of the system can be decomposed into three parts, 

+ AQ kin . (1) 

In the thin-wall approximation the bulk and surface terms are expressed through the bubble 
radius R in a standard way, 

Afi^ fc = --R 3 (P H - Pq) , A{l surf = 4vrP 2 7 , (2) 
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where Pr and Pq are the pressures of the bulk H and Q phases, 7 is the effective surface 



where A£ = Sq — £h is the difference of energy (more exactly, enthalpy) densities of the 
two bulk phases. It is evident that this term is negative, since the dense Q phase is replaced 
in the bubble by the dilute H phase (typically, £q ^> £u). Thus, the bubble formation is 
favoured by the collective expansion. Moreover, the nucleation can start now even above 
the critical line, when P H < Pq, and the standard theory would predict no growing bubbles. 
In principle, the phase separation can start as early as the metastable H state appears in 
the thermodynamic potential, and a stable interface between the phases may exist. 

Using Eqs. (0) and (|3|) one can determine the critical bubble radius R c , corresponding 
to the top of the potential barrier in AQ(R). The condition d^AQ = leads to a cubic 
equation for R c . When H — > the kinetic term vanishes and this equation gives a standard 
Laplace formula for the critical bubble ||. However for realistic parameters (see below) the 
kinetic term dominates. In particular, in the vicinity of the critical line, when Ph ~ Pq, one 
can consider the bulk term perturbatively. Then one obtains 



The bubbles with R > R c will expand further while those with R < R c will eventually shrink. 
From the above consideration one should conclude that in a rapidly expanding system an 
appreciable amount of nucleation bubbles and even empty cavities will be created already 
above the critical line. 

The bubble formation and growth will also continue below the critical line. Previously 
formed bubbles will now grow faster due to increasing pressure difference, Ph — Pq > 0, 
between the two phases. It is most likely that the conversion of Q matter on the bubble 
boundary is not fast enough to saturate the H phase. Therefore, a fast expansion may lead 



tension. The last term in Eq. ([!]) accounts for the change in the local kinetic energy of 



expanding matter. In the same approximation it can be evaluated as 




(3) 
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to a deeper cooling of the H phase inside the bubbles compared to the surrounding Q matter. 
Strictly speaking, such a system cannot be characterized by the unique temperature. At 
some stage H bubbles will percolate, and the topology of the system will change. Now 
isolated regions of the Q phase (Q droplets) will be surrounded by the undersaturated vapor 
of the H phase. 

Standard thermodynamical concepts cannot be used in this non-equilibrium situation. 
However, the characteristic droplet size can be estimated by applying the energy balance 
consideration first proposed by Grady [pT0| , pLl[1 in the study of dynamical fragmentation. 
The idea is that the fragmentation of expanding matter is a local process minimizing the 
sum of surface and kinetic (dilational) energies per fragment volume. The predictions of 
this simple model are in reasonable agreement with molecular dynamics simulations |TT|JT^] 
and experimental data on dynamical fragmentation of fluids and solids (see e.g. [ p~0| , p~T| ,p~3f| ) . 
As shown in Ref. [f[4 |, this prescription works fairly well also for multifragmentation of 



expanding nuclei, where the standard statistical approach fails. 

Let us imagine an expanding spherical Q droplet embedded in the background of the 
dilute H phase. In the droplet rest frame the change of thermodynamic potential compared 
to the uniform H phase is given by the same expression ([I]) but with indexes H and Q 
interchanged. The kinetic term is positive now. According to the Grady's prescription, the 
characteristic droplet radius, R*, can be determined by minimizing 

(^) = - OPq - Pn) + ^ + ^SH^ . (5) 

\ V J droplet ri 1U 

It is worth noting that the collective kinetic energy term acts here as an effective long-range 
potential, similar to the Coulomb potential in nuclei. Since the bulk term does not depend on 
R the minimization condition constitutes the balance between the collective kinetic energy 
and interface energy. This leads to an optimum droplet radius 



\ASH 2 , 

It should be noticed that this radius is expressed through the same combination of model 
parameters as the critical bubble radius at Ph ~ Pq, Eq- @, but with a slightly bigger 



numerical coefficient. This suggests that in vicinity of the critical line the H and Q phases 
occupy roughly equal fractions of the total volume. Fast expansion (large H) may lead to 
very small droplets. This mixed state of matter is far from thermodynamical equilibrium 
because of the excessive interfacial energy and undersaturation of the H phase. One can 
say that the metastable Q matter is torn apart by a mechanical strain associated with the 
collective expansion. This has a direct analogy with the fragmentation of pressurized fluids 
leaving nozzles [OjIH]- m a similar way, splashed water forms droplets which have nothing 
to do with the equilibrium liquid-gas phase transition. 

At ultrarelativistic collision energies associated with RHIC and LHC experiments, the 
expansion will be very anisotropic, with its strongest component along the beam direction. 
Applying the same consideration for the anisotropic flow, one can see that the characteristic 
size of inhomogeneities in each direction is determined by the respective Hubble constant 
acting in this direction. Thus, the resulting structures will have smaller size in the direction 
of stronger flow. Therefore, in the case of strong one-dimensional expansion the inhomo- 
geneities associated with the phase separation will rearrange themselves into pancake-like 
slabs of Q matter layered by the dilute H phase. The characteristic width of the slab is given 
by Eq. (||) with a slightly different geometrical factor. At a later stage the slabs will further 
fragment into smaller droplets. In general one should expect the emergence of complicated 
multi-fractal structures. 

The driving force for expansion is the pressure gradient, VP = ci?V£, which depends on 
the sound velocity in the matter, c s . In the vicinity of the critical line one may expect a 
"soft point" |T5| , p!6| where c s is smallest and the ability of matter to generate the collective 
expansion is minimal (small H). If the initial state of the Q phase is close to this point, 
the primordial bubbles or droplets will be biggest. Increasing initial pressure will result in 
a faster expansion and smaller droplets. For numerical estimates I choose two values of the 
Hubble constant: if~ 1 =20 fm/c to represent the slow expansion from the soft point ]T3[ and 
H~ 1 =6 fm/c for the fast expansion ||. 

One should also specify two other parameters, 7 and AS. The surface tension 7 is a 
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subject of debate at present. Lattice simulations indicate that at the critical point it could 
be as small as a few MeV/fm 2 . However, for the non-equilibrium situation discussed here 
the values of 10-20 MeV/fm 2 , which follow from effective chiral models, should be more 
appropriate. As a compromise, the value 7 = 10 MeV/fm 2 is used below. It is clear that 
AS should be close to the latent heat of the transition, i.e. about 0.5 -j- 1 GeV/fm 3 . One 
can also estimate AS by realizing that nucleons and heavy mesons are the smallest droplets 
of the Q phase. For estimates I take AS = 0.5 GeV/fm 3 , i.e. the energy density inside the 
nucleon. Substituting these values in Eq. @ one gets R*=3A fm for H~ 1 =20 fm/c and 
R*=1.5 fm for H~ x =^ fm/c. 

In the lowest-order approximation the characteristic droplet mass is M* pa ASV . For 
spherical and slab-like droplets one gets respectively 

JC»lTj£i. ^-2S(A£)^(^) 1/3 , (7) 

where S is the slab transverse area. It is interesting to note that in this approximation M* p 
is independent of AS. For the two values of R* given above M* p is ~100 GeV and ~10 
GeV, respectively. The slab-like structures would have even larger mass, since S could be 
of order of the transverse system size. Using the minimum information principle one can 



show |TT| , ^4| that the distribution of droplets should follow an exponential law, exp ( — 



M 
M* 



Therefore, with 1% probability one can find droplets as heavy as 5M*. 

After separation, the droplets recede from each other according to the global Hubble 
expansion, predominantly along the beam direction. Therefore, their center-of-mass rapidi- 
ties are in one-to-one correspondence with their spatial positions. Presumably they will be 
distributed more or less evenly between the target and projectile rapidities. At this late 
stage it is unlikely that the thermodynamical equilibrium will be re-established between the 
Q and H phases or within the H phase alone. If this were to happen, the final H phase would 
be uniform, and thus there would be no traces of the mixed phase in the final state. 

The final fate of individual droplets depends on their sizes and on details of the equation 
of state. Due to the counter-acting pressure of the H phase and additional Laplace pressure 



the initial expansion will slow down. In smaller droplets the expansion and cooling may even 
reverse to the contraction and reheating. The conversion of Q droplets into the H phase 
may proceed through formation of a deflagration front |16| or evaporation of hadrons from 
the surface ]T7|] . Bigger droplets may expand further until they enter the region of spinodal 
instability. As shown in Ref. [|18|| , the characteristic time of the "rolling down" process is 
rather short, ~ 1 fm/c, so that the Q droplets will be converted rapidly into the H phase. 
The energy released in this process can be transferred partly into the collective oscillations 
of the (a, 7r) fields. Numerical simulations Jl9|,^0j show that these oscillations persist for 



a long time and give rise to soft pion radiation. One should also expect the formation of 
Disoriented Chiral Condensates (DCC) in the voids between droplets. 

Since rescatterings in the dilute H phase are rare, most hadrons produced from individual 
droplets will go directly into detectors. One may guess that the number of produced hadrons 
is proportional to the droplet mass. Each droplet will give a bump in the hadron rapidity 
distribution around its center-of-mass rapidity |L8|]. If emitted particles have a Boltzmann 



spectrum, the width of the bump will be Sy ~ 2JT/m, where T is the droplet temperature 
and m is the particle mass. At T ~ 100 MeV this gives Sy ~ 2 for pions and 5y ~ 1 for 
nucleons. These spectra might be slightly modified by the residual expansion of droplets 
and their transverse motion. The resulting rapidity distribution in a single event will be a 
superposition of contributions from different droplets, and therefore it will exhibit strong 
non-statistical fluctuations. The fluctuations will be more pronounced if primordial droplets 
are big, as expected in the vicinity of the soft point. If droplets as heavy as 100 GeV 
are formed, each of them will produce up to ~300 pions within a narrow rapidity interval, 
Sy ~ 1. Such bumps can be easily resolved and analyzed. Critical fluctuations of similar 
nature were discussed recently in Ref. [[H] . 



Some unusual events produced by high-energy cosmic nuclei have been already seen 
by the JACEE collaboration P^[ . Unfortunately, they are very few and it is difficult to 
draw definite conclusions by analyzing them. We should be prepared to see thousands 
of such events in the future RHIC and LHC experiments. It is clear that the nontrivial 
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structure of the hadronic spectra will be washed out to a great extent when averaging 
over many events. Therefore, more sophisticated methods of the event sample analysis 
should be used. The simplest one is to search for non-statistical fluctuations in the hadron 
multiplicity distributions measured in a fixed rapidity bin ||23|| . One can also study the 
correlation of multiplicities in neighbouring rapidity bins, bump-bump correlations etc. Such 



standard methods as intermittency and commulant moments wavelet transforms |24 



HBT interferometry |25j can also be useful. All these studies should be done at different 



collision energies to identify the phase transition threshold. The predicted dependence on 
the Hubble constant and the geometry of reaction, Eq. (0), can be checked in collisions with 
different ion masses and impact parameters. 

One should bear in mind two important points. First, if the expansion trajectory goes 
close to the (tri) critical point, both 7 and 5£ will tend to zero and the critical fluctuations 
will be less pronounced. Second, if a first order phase transition is possible only in the 
baryon-rich matter, then the Q droplets should have much higher baryon density than the 
hadronic phase ]3[]. In this case one should expect strong non-statistical fluctuations in the 
distribution of the net baryon charge. 

In conclusion, simple arguments based on the homogeneous nucleation picture and the 
energy balance consideration demonstrate that a first order phase transition in rapidly ex- 
panding matter should proceed through a non-equilibrium stage when the metastable phase 
fragments into droplets. If QCD matter undergoes a first order phase transition, it will 
manifest itself in relativistic heavy-ion collisions by the formation of droplets of quark-gluon 
plasma. The primordial droplets should be biggest in the vicinity of the soft point where the 
expansion is slowest. The fragmentation of plasma might be accompanied by the formation 
of multiple DCC domains and enhanced soft-pion radiation. Subsequent hadronization of 
QGP droplets will lead to large non-statistical fluctuations in the hadron rapidity spectra 
in individual events. All these novel phenomena can only be detected through dedicated 
event-by-event analysis of experimental data. 
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FIG. 1. Schematic view of the effective thermodynamic potential per volume fl/V as a function 
of the order parameter field a at ir = 0, as predicted by the linear u-model in the chiral limit 
= [1]. The curves from bottom to top correspond to the different stages of the isentropic 
expansion of homogeneous matter starting from T=100 MeV and ^=750 MeV (curve 1). The 
upper curve 5 is the vacuum potential. The other curves are discussed in the text. 
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